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THEORETICAL PROPERTIES OF THE WEIGHTED FELLER-PARETO 
AND RELATED DISTRIBUTIONS 


OLUSEYI ODUBOTE AND BRODERICK O. OLUYEDE 


ABSTRACT. In this paper, for the first time, a new six-parameter class of distributions called 
weighted Feller-Pareto (WFP) and related family of distributions is proposed. This new class 
of distributions contains several other Pareto-type distributions such as length-biased (LB) 
Pareto, weighted Pareto (WP I, IJ, III, and IV), and Pareto (P I, II, III, and IV) distributions 
as special cases. The pdf, cdf, hazard and reverse hazard functions, monotonicity properties, 
moments, entropy measures including Renyi, Shannon and s-entropies are derived. 


1. INTRODUCTION 


Pareto distributions provide models for many applications in social, natural and physical 
sciences and are related to many other families of distribution. A hierarchy of the Pareto 
distributions has been established starting from the classical Pareto (I) to Pareto (IV) distri- 
butions with subsequent additional parameters related to location, shape and inequality. A 
general version of this family of distributions is called the Pareto (IV) distribution. Pareto 
distribution has applications in a wide variety of settings including clusters of Bose-Einstein 
condensate near absolute zero, file size distribution of internet traffic that uses the TCP 
protocol, values of oil reserves in oil fields, standardized price returns on individual stocks, 
to mention a few areas. 

Brazauskas (2002) determined the exact form of the Fisher information matrix for the 
Feller-Pareto distribution. Rizzo (2009) developed a new approach to goodness of fit test 
for Pareto distributions. Riabi et al. (2010) obtained entropy measures for the family of 
weighted Pareto-type distributions with the the general weight function w(x; k,t,i, j) = 
xte F(x) FI (x), where F(x) and F(x) = 1— F(a) are the cumulative distribution function 
(cdf) and survival or reliability function, respectively. 


1.1. Some Basic Utility Notions. Suppose the distribution of a continuous random vari- 
able X has the parameter set 6* = {01, 62,--- ,@,}. Let the probability density function (pdf) 
of X be given by f(x; 6*). The cumulative distribution function (cdf) of X, is defined to be 
F(x;0*) = f7 f(t;0*) dt. The hazard rate and reverse hazard rate functions are given by 





h(x; 0*) = H and T(x; 0*) = fen, respectively, where F'(x; 6*) is the survival or reli- 
ability function. The following useful functions are applied in subsequent sections. The gam- 
ma function is given by P(x) = fy t?~'e~ dt. The first and the second derivative of the gam- 


ma function are given by: I(x) = fy tt (logt)e™ dt, and I(x) = fy t®-*(log t)?e™ dt, 
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respectively. The digamma function is defined by U(x) = I’(x)/I'(x). The lower incom- 
plete gamma function and upper incomplete gamma function are y(s, x) = ie ts-te—t dt and 
I'(s,z) = f° tte dt, respectively. 
Let a,b > 0, then 

P (a)T (b) 


(1) i. z**(1— 2)! In(x) dz = Teci 


1.2. Introduction to Weighted Distributions. Statistical applications of weighted dis- 
tributions, especially to the analysis of data relating to human population and ecology can 
be found in Patil and Rao (1978). To introduce the concept of a weighted distribution, sup- 
pose X is a non-negative random variable (rv) with its natural probability density function 
(pdf) f(x;0), where the natural parameter is 0 € Q (Q is the parameter space). Suppose a 
realization x of X under f(x; 6) enters the investigator’s record with probability proportional 
to w(x; 2), so that the recording (weight) function w(x; 8) is a non-negative function with 
the parameter 8 representing the recording (sighting) mechanism. Clearly, the recorded x is 
not an observation on X, but on the rv Xu, having a pdf 

(2) fola 0,p) = MPN) 

where w is the normalizing factor obtained to make the total probability equal to unity 
by choosing 0 < w = Elw(X,()| < oo. The random variable X,, is called the weighted 
version of X, and its distribution is related to that of X. The distribution of X, is called the 
weighted distribution with weight function w. Note that the weight function w(x, 3) need not 
lie between zero and one, and actually may exceed unity. For example, when w(x; 3) = x, in 
which case X* = Xw is called the size-biased version of X. The distribution of X* is called 
the size-biased distribution with pdf f*(x;0) = aoa where 0 < u = E|X] < œ. The pdf 
f* is called the length-biased or size-biased version of f, and the corresponding observational 
mechanism is called length-biased or size-biased sampling. Weighted distributions have seen 
much use as a tool in the selection of appropriate models for observed data drawn without a 
proper frame. In many situations the model given above is appropriate, and the statistical 
problems that arise are the determination of a suitable weight function, w(x; 3), and drawing 
inferences on 6. Appropriate statistical modeling helps accomplish unbiased inference in spite 
of the biased data and, at times, even provides a more informative and economic setup. See 
Rao (1965), Patel and Rao (1978), Oluyede (1999), Nanda and Jain (1999), Gupta and 
Keating (1985) and references therein for a comprehensive review and additional details on 
weighted distributions. 

Motivated by various applications of Pareto distribution in several areas including reliabil- 
ity, exponential tilting (weighting) in finance and actuarial sciences, as well as in economics, 
we construct and present some statistical properties of a new class of generalized Pareto-type 
distribution called the Weighted Feller-Pareto (WFP) distribution. 

The aim of this paper is to propose and study a generalization of the Pareto distribution 
via the weighted Feller-Pareto distribution, and obtain a larger class of flexible parametric 
models with applications in reliability, actuarial science, economics, finance and telecommu- 
nications. This paper is organized as follows. Section 2 contains some utility notions and 
basic results. The weighted Feller-Pareto distribution is introduced in section 2, including 
the cumulative distribution function (cdf), pdf, hazard and reverse hazard functions and 
monotonicity properties. In section 3, moments of the WFP distribution are presented. The 


(U(a) — U(a +). 
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mean, variance, standard deviation, coefficients of variation, skewness, and kurtosis are read- 
ily obtained from the moments. Section 4 contains measures of uncertainty including Renyi, 
Shannon and s—entropies of the distribution. Some concluding remarks are given in section 
5. 


2. THE WEIGHTED FELLER-PARETO CLASS OF DISTRIBUTIONS 


In this section, the weighted Feller-Pareto class of distributions is presented. First, we 
discuss the Feller-Pareto distribution, its properties and some sub-models. Some sub-models 
of the FP distribution are given in Table 1 below. 


2.1. Feller-Pareto Distribution. In this section, we take a close look at a more general 
form of the Pareto distribution called the Feller-Pareto distribution which traces it’s root 
back to Feller (1971). 


Definition 2.1. The Feller-Pareto distribution called FP distribution for short, is defined 
as the distribution of the random variable Y = u + 0(X™! — 1V, where X follows a beta 
distribution with parameters a and 8, 0 > 0 and y > 0, that is, Y ~ FP(u,0, Y,a, B) if the 
pdf of Y is of the form 


D teina = aal EAT 





Bla, B)Oy\ 8 0 


for —o0 < u < 00, a > 0, B > 0, 0 > 0, y > 0, and y > p, where B(a, 8) = WER. 
The transformed beta (TB) distribution is a special case of the FP distribution. The pdf 


of the transformed beta distribution is given by 
1 78 
(2/0) yi 
Bia, B)ax [1 + (x/0)7]e+8 
Therefore, from (4), TB(0,y,a, 8) = FP(0,0,1/7, a, 8). The family of Pareto distributions 


(Pareto I to Pareto IV) can be readily obtained for specified values of the parameters pu, 0, 
y, a, and $. 





(4) frep(z; O50) = 


TABLE 1. Some Sub-Models of the FP Distribution 








Family name Symbol Density function 


= (an) 
FPO) | FPCyin.81,0,1)| a a + ($4) | 


I —(a+1) 
FP(I) FP(y; p, 0, l,a, 1) A (5) f + (=) | 

















FPI | FP(y;#,6,7,1,1) F 




















T=] To —(a+l) 
FP(IV) | FP(y;11,,7,051) ml) hefy 
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2.2. Moments of Feller-Pareto Distribution. The k moment of the random variable 
YE under FP distribution is given by: 


yay k 7 œ 1 ea k+8—1 msi 44 -(a+8) 
( 7 ) 7 | Fa a 7 ) + ( 9 ) | L 


1 


y ae 
Let (=) =z 40<t<i1, then ty = b( 34 :) aa and 


r y — u k 7 [ by t y(k+£-1) 1 —(a+8) t y= dt 
0 Jy Bla, B)Oy\1—-t 1-¢t 1-t (1 — t)? 


= = i ky+B-1 _ 4\a-ky-1 
= a (1 -t) dt 
Pky + L) (a — ky) 
(a) P'() l 


for k = 0,1,...., a — ky # 0, —1,—2,...., and a — ky > 0. If u = 0, then the kt? moment 
reduces to 














O"T (ky + B)P (a — ky) 
T'(a)P(B) 


for = <k< F Note that if u = 0 and 0 = 1, then we have, 





E(Y*) = 


P(ky + B)E(a — ky) 
P(a)P(S) 
k = 0,1,2,..., a — ky 4 0,-1,—-2,..... The mean, variance, coefficient of variation (cv), 


coefficient of skewness (cs) and coefficient of kurtosis (ck) of the Feller-Pareto distribution 
(when u = 0) are given by 


Aly + Ba — ¥) 


BY) = 





























ce oC) (mae 
52, — PERT + BT (a = 27) _ as + AT (a — a 
i F(a) (8) r(a)r(8) ' 
Pr Crer (a27) — (Soare) 
ay vei TrA) 
Be — 7) 
Tari) 
E(Y3) — 340? — p? 
CSpp = r 5 3H a 
E(Y*) — 44 E(Y) + 6P E(Y?) — 4u’ E(Y) + ut 
CK rp = ot ? 
where u = upp, o = \/o%p, E(Y’) = ign and E(Y*) = HESE, See Arnold 


(1983), Luceno (2006) and references therein. 
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2.3. Weighted Feller-Pareto Distribution. In this section, we present the weighted 
Feller-Pareto (WFP) distribution. Some WFP sub-models are also presented in this sec- 
tion. Table 2 contains the pdfs of the sub-models for the WFP distribution. First, consider 
the weight function w(y;k) = y*. The WFP pdf fwrp(y), when u = 0 and 6 = 1, is given 


by 
wen) 2 ~ OT re — ky) one f i (5) i 


for = < k < ©. The length-biased (k = 1) Feller-Pareto (LBFP) pdf is 


aroi (8) OT 


k 
for y < a and u = 0. In general, with the weight function w(y; u, 0, k) = (=) , we obtain 
the WFP pdf as: 














(6) fwer(y) = ee a h+ (9 
D(ky + B)P(a— ky)B(a, B)by\ 0 0 
for k = 0,1,2,...; y > u, & — ky # 0,—1,—2,....... , and a — ky > 0. The cdf of WFP 
distribution is given by: 
k+8-—1 
: ['(a + p) 2 
F(y; a, 8, htk) i TE 


d 
+ky)T (a — ky) Farg 
pe (=) 
(44#)G/y) 
Ta tt B)B (Ee ky f p, a + y= ky) 
6 — 
j +A h) 
for a—ky £ 0,-1, -2,......,a+y—ky 4 0, -1, -2,......, where B(t; a,b) = i. y? 1(1-y)>ldy 
is the incomplete beta function. The plots of the pdf and cdf of the WFP distribution in 


Figures 1 and 2 suggest that the additional parameters k controls the shape and tail weight 
of the distribution. 








Weighted Feller-Pareto Pdf Weighted Feller-Pareto Pdf with K=1 

















—— PI 0.5 6-0.5 0.2 a=0.4 B=1 
Pll w=0.6 60.4 y=0.2 a=0.4 B=1 
—— PII p=0.7 60.5 70.2 0=1B=1 H 
—— PIV p=0.7 60.6 0.5 a=0.8 B=1 


10.8 60.5 Y=0.2 o=0.9 B=O.4 k=0 
—— 170.8 60.5 0.2 o=0.9 B=O.4 k=1 
—— 70.8 00.5 y=0.2 00.9 B=O.4 k=2 j 0.025 
0.8 60.5 90.2 o=0.9 B=O.4 k=3 
—— u08 60.5 90.2 o=0.9 B=O.4 k=4 
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FIGURE 1. PDF of the Weighted Feller-Pareto with k=1 and different values of k 


FIGURE 2. CDF of the Weighted Feller-Pareto with k=1 and different values of k 


CDF WFP 
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o=0.8; B=0.6; y=0.2; p=0; 60.3; k=0; 
—:—'-g=0.8; B=0.6; y=0.2; p=0; 60.3; k=1; 
o=0.8; B=0.65; y=0.2; p=0; 6=0.3; k=2; 


+ 

















TABLE 2. Sub-Models of the WFP Distribution 








Family 
name 


Symbol 


Density function 





WP(I) 


FP(6,0,1,0,1,k 


Tr(a) 
T(ky+1)0(a—k)B(a,1)0 





h 





WP(II) 


FP(u,0,1,a,,k) 





1 





T(a) l 
T(ky+1)0 (a—k) B(a,1)0 





WP(III) 


P(e 07,4, lk 


T 
y 
YK 


15-2 
T(ky+1)lG—ky) B11) 6y : 0 








WPUIV) 





FP(p, 0, Y, Q, 1, k 





yY 





F 
a= T(a+1) 





T(ky+1)l(a-ky)0y 





I 
1 | 





2.4. 


given by: 


hpr(y; k, a, 8, u, 7,9, k) 


Hazard and Reverse Hazard Functions. In this section, hazard and reverse hazard 
functions of the WFP distribution are presented. Graphs of the hazard function for selected 
values of the model parameters are also given. The hazard and reverse hazard functions are 














f(y; a, b, p, 7,9, k) 
1— F(y; a, B, u, V0 k) 
k+-1 al kom 
sift (St) l) 
O40 (6+ky)T(a—ky) 
T(a+8) (G4)G/)) 
l= T+A a- 8 (a +E amik + ba +y- ky 
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and 


f(y; a, Bb, 7, 9, k) 
F(y; a, By, 7,0,%) 
By 1 


EE k+5 =) 22 
O90 (B+ky) r (a—ky) (=) h T (#) | 


T'(a+8) (13E) g , 
Se P (margi hy + Ba + k 


for a— ky # 0, —1, —2,...... a+ y—ky £0,-1,-2,....... , respectively. Graphs of the WFP 
hazard rate function for selected values of the model parameters are given in Figure 3. The 
graphs shows unimodal and upside down bathtub shapes. 





Tr(y; a, b, H, V5.0, k) = 














—:—:-a=0.8; B=0.6; y=0.2; p=0; 60.3; k=0; 

o=0.8; B=0.6; y=0.3; p=0; 60.3; k=1; 

~~~ o=0.9; B=0.6; y=0.4; p=0; 0.3; k=2; 
s 








FIGURE 3. Hazard function of the Weighted Feller-Pareto with k=0 and dif- 
ferent values of k 


2.5. Monotonicity Properties. In this section, monotonicity properties of the WFP dis- 
tribution are presented. The log of the WFP pdf is given by: 


In(fwree(y)) = nlogr(a+ 8) —nlogI(ky + p)—nlogr(a — ky) — nlog 6y 
+ lypsi log =E — (a + 8) log |1 + EEN 
y 0 0 
Now differentiating In(fwrp(y)) with respect to y, we have 


Olnfwrrly) _ (k -1+ 2) E (a + 8) (2) g 


Oy y— UL O N7 
vol (=) 
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and solving for y, we have 


Olnfwrp(y;9 VQ ay p, os y=0(7 =p = oya 














Oy KY = y-a 
Note that 
Olnfwrr(y;9 V Q ay p, k) Y= Pky i 
Dy <0 <4 y >90 Prensa T H, 
and 








ln fwrply; 0 V QA a, p, 2 >o y <e (1E) 
Oy 


5 
The mode of the WFP distribution is given by y, = 0 (g) + pL. 


3. MOMENTS OF WFP DISTRIBUTION 


Recall that the kt? moment of the random variable XE under the FP distribution are 
given by 





n” - 2) _ Pky + A)T (a - ky) 

0 F(a) l 
for k = 0,1,...., a— ky #0, —1, —2,....,a — ky > 0. Now, we derive the rt” moments of the 
random variable Xo under WFP distribution. This is given by 





Om — T 


y-1 
tet (1) = 1h. 0< £< 1, then dy = 07( 35 ) Tope and 


a Í J j e een (i L) per 


(a+) y-1 
t t 1 
P+] Ga) eal 


2 r(a)Tr(8) [ {eR EC = mea ae ar 
P(ky + B)P(a — ky) Bla, p) Jo 
Dky+B)P(a-ky) 
for a — ky # 0,-1, -2,...... ,a—ry # 0,-1,-2,...... „œa — y(r + k) # 0,-1,-2,...... eS 
0,1,2,..... The mean, variance, coefficients of variation (cv), skewness (cs) and kurtosis (ck) 
can be readily obtained from equation (7). 








x 
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4. MEASURES OF UNCERTAINTY FOR THE WEIGHTED FELLER-PARETO DISTRIBUTION 


The concept of entropy was introduced by Shannon (1948) in the nineteenth century. 
During the last couple of decades a number of research papers have extended Shannon’s 
original work. Among them are Park (1995), Renyi (1961) who developed a one-parameter 
extension of Shannon entropy. Wong and Chen (1990) provided some results on Shannon 
entropy for order statistics. The concept of entropy plays a vital role in information theory. 
The entropy of a random variable is defined in terms of its probability distribution and can 
be shown to be a good measure of randomness or uncertainty. In this section, we present 
Renyi entropy, Shannon entropy and s—entropy for the WFP distribution. 


4.1. Shannon Entropy. Shannon entropy (1948) for a continuous random variable Y with 
WFP pdf fwrp(y) is defined as 





(8) E(-los(fweetY)) = f “(Besa ie. 
Note that 
E ey) -ios ( ern | | log? (a) + log? (B) — logl (ky + 8) 


logl (a — ky) + (: + i — 1)roa( 27+) 


= (a+ Bjiog|a + (232) ]. 


Now, Shannon entropy for the weighted Feller-Pareto distribution is 





setin = -ioa Sr | + log (a) + logl (8) — log (ky + 2) 
~ satiety) a fo(“S2)) 
(9) - (a+ 96 [toaf1+ (ZEY. 


xy 
Now, with the substitution (=) = ;4 for 0 < t < 1, we can readily obtain both 


E toy (s)| and E toy j+ (=) | | , so that Shannon entropy for the WFP distribution 
is given by 
E(-logfwrp(Y)) = —logB(a, p) + log(@y) + P(a) + (8) — log (ky + 8) 
log’ (a — ky) + (ky + 8 — ¥)[¥(8) — y(a)] 
(a + B)((a) — Ya + §)), 


(10) 


where Ņ(.) = TE is the digamma function. 


10 ODUBOTE AND OLUYEDE 


4.2. s-Entropy for Weighted Feller-Pareto Distribution. The s-entropy of the WFP 
distribution is defined as 
1 CO ps ; 
H,(fwrp) =% h = iB fives) ifs #1, s> 0, 
E(-logfwrp(Y)), ifs=1. 
Note that 


k mE I In i aie “pIE By) ( z T 





y-1 
. . . = d 
and using the substitution (=) = for 0 <t < 1, so that dy = n1) ao? we 


get 
CO 1 1 S 
g dy = i ponent” =e a OTH — t) tdt 
J, edis = | agnt oe] ea- 


all (a + p) (ks + s8 — sy +y) (sa — ks + sy- y) 


(07) [[(ky + (a — ky)sT'(s8 + sa) 





(11) 
for s > 0, s # 1. Consequently, s-entropy for the WFP distribution is given by 
1 _, (at b) Tks + s8 — sy +y)I(sa — ks + sy-— 
H.(fwer) = 5l _ (oy ECAT (ks + s8 — oy + VE daa 
S (ky + Ara — ky) PT (SB + sa) 
for S > 0, S + 1, a—ky # 0, =], = 2, e.. sa+sy—ks—y Æ 0, =I, —2, sees ks+sB—sy+y > 0, 
and sa — ks + sy—7 > 0. 





4.3. Renyi Entropy. Renyi entropy (1961) for the WFP distribution is presented in this 
section. Note that Renyi entropy is given by 


(12) Hy(Fvee) = pag ( f Uner(e'de), 9> 0,941 


From equation (10), we obtain Renyi entropy as follows: 
fs nee 1 log [Es + A) T(ks + 88 — sy +7) (sa — ks + sy — 7) 
A (07) (ky + BI (a — ky)]'T (s8 + sa) 
for s > 0, s # 1, a— ky # 0, —1, —2, ..., sa+sy— ks—y + 0, —1, —2, ...., ks tsb — sy +y > 0, 
and sa — ks + sy — y > 0. 











5. CONCLUDING REMARKS 


In this paper, a new six-parameter class of distributions called weighted Feller-Pareto 
(WFP) distribution is constructed and studied. The pdf, cdf, hazard and reverse hazard 
functions, monotonicity properties are presented. Measures of uncertainty including Renyi, 
Shannon and s-entropies are derived. 
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TABLE 3. Shannon and s-entropies of the sub-models of the FP distribution 

















Family Shannon entropy s-entropy 
name 
Tat T(a+1)]§T(stas— 
P(I) log ( st) P(a)-1+(a+1)(W(a)— | 2 * J1 — oE T etas) 
Yla +1)) i : 
1)]'T(s+as—7) 














P(IL) log( "32 Tia) —1 + (a +| J1- oE 


















































1)y(a) = yla + 1) 
P(IIT) lg( 52) 24 2(4(1) y(2)) xli e 
P(IV) log( Se n) P(a)-1-(1-7)[v(1)- | 4* |1 fy Neen 
¥(a)] + (a+ D(a) - yla + 1)) ; 





TABLE 4. Shannon and s-entropies of the WFP Distribution 








Family name | Shannon entropy s-entropy 


WP(D) logB(a, 1) +log(8)—T(a)—1+logl (k "E tee] 
1)+logr (a — k) — kiy (1)—y(a)]+ (a 

1)(W(a) — yla + 1)) 

W P(II) logB(a, 1) +log(0)—T'(a)—1+logl (k Ah Toa 
1)+logr (a — k) —k[(1) —v(a)] + (a 

Hlo) -ya D) 

WP(III) logB(1, 1) +log(67) —2+logl(ky+1)+ | 2 |: RER 
logl (1—ky)—(ky+1—-7)[¥()—-v(1)] 

2(y(1) — ¥(2)) 

WPUV) logB(a,1) + log(@y) — F(a) — 1 Ah T 

log! (ky+1)+logT(a— ky) — ( 


J \ | / | [F (a+1)] Tr (ksy+y) (sa—ksy+s—y) 
NIPDAU) | eenas] 
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